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祁连山 及 周边 降水 分 布 聚 类 检验 和 典型 流域 增 雨 效果 评价 
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摘 要 : 祁连山 及 周边 区 域 地 形 复杂 ,降水 分 布 极 不 均匀 ,如 何 科学 地 分 区 分 型 ,对 把 握 研究 区 域 降 
水 分 布 特征 具有 极为 重要 的 意义 。 利 用 该 区 域 31 个 国家 基本 站 1961—2020 年 $ 一 9 月 降水 量 资 
料 , 采 用 主 成 分 分 析 方 法 (PCA) 对 该 区 域 降 水 量 进 行 分 析 , 再 通过 聚 类 分 析 (CAST) 对 该 结果 进行 


显著 性 检验 ,最 后 将 分 区 分 型 结果 应 用 于 人 工 增 月 


作业 效果 评估 。 结 果 表 明 :(1) 祁连山 及 周边 可 


分 成 7 个 区 域 (Z1~Z7) ,累积 方差 贡献 率 超过 78%。(2) 这 7 个 区 域 分 别 以 乐 都 、 海 曼 、 野 牛 沟 .武威 、 
高 合 、 临 泽 、 刚 察 为 中 心 点 ,划分 结果 与 降水 量 、 地 理 地 形 、 海 拔高 度 显著 关联 ,(3) 依托 区 域 历史 回 
归 统 计 方 法 对 人 工 增 雨 作业 效果 进行 评价 ,发 现 1992 一 2020 年 石 羊 河流 域 5 一 9 月 绝对 增 雨 量 、 相 
对 增 雨 率 分 别 为 8.91 mm、6.51%, 其 中 7 月 最 高 (6.30 mm、21.86%),8 月 次 之 (5.44 mm、16.11% ) 。 


基于 地 


降水 量 的 作业 效果 评价 往往 受 对 比 区 选择 的 影响 , 聚 类 检验 方法 不 仅 有 助 于 复杂 地 形 下 


降水 量 的 分 区 研究 ,还 为 科学 选择 对 比 区 提供 了 客观 依据 。 
X 键 词 : 聚 类 检验 ; 主 成 分 分 析 ; 分 区 分 型 ; 增 雨 效果 评价 


文章 编号 : 


祁连山 地 处 青藏 . 蒙 新 .黄土 高 原 交 汇 地 人 带 ,位 
居 我 国 干旱 半 和 干旱 区 内 陆 腹 地 , 横 跨 青海 省 东北 部 
与 甘肃 省 西部 地 区 。 整 个 祁连山 区 属 裙 皱 断 块 , 自 
西北 向 东南 依次 为 大 雪山 、 托 来 山 、 托 来 南山 、 野 马 
南山 AA L 党 河南 山 . 土 尔 根 达 坂 山 、 柴 达 木 
山 和 宗 务 隆 山 ,平均 海拔 4000 m 以 上 ,山脉 东西 长 
1000 km ,南北 宽 300 km ,范围 广大 ,地 形 复杂 ,降水 
区 域 性 分 布 特征 尤为 明显 。 作 为 西北 干旱 区 的 一 
条 湿 舌 ,祁连山 不 仅 是 诸多 河流 与 水 系 的 发 源 地 ， 
也 是 周边 区 域 生存 与 发 展 的 屏障 。 祁 连 山北 侧 是 
石 羊 河 ` 黑 河 ` 玻 勒 河 三 大 内 陆 河 径流 区 ,是 河西 绿 
洲 存活 的 基础 ;祁连山 南 侧 是 大 通 河 、 涅 水 河流 经 
地 ,是 河 涅 流域 生存 的 命脉 。 祁 连 出 的 存在 影响 着 
周边 大 范围 地 区 ,形成 一 个 相互 影响 且 不 可 分 割 的 
整体 。 近 年 来 祁连山 雪线 上 升 、 冰 川 面积 锐 减 .水 
资源 严重 短缺 等 问题 日 益 凸 显 ,成 为 社会 经 济 与 生 


出 现 , 探 讨 祁连山 及 周边 降水 分 布 与 变化 规律 有 着 
重要 的 现实 意义 。 

人 工 增 雨 是 开发 祁连山 空中 云 水 资源 以 及 组 
解 水 资源 短缺 的 有 效 手 段 之 一 ,其 有 效 性 与 科学 性 
已 经 被 越 来 越 多 的 研究 者 所 重视 ”。 作 为 人 工 增 
雨 试验 研究 的 技术 难点 ,作业 效果 评价 是 一 项 极其 
复杂 而 重要 的 工作 。 由 于 自然 降水 变 率 极 大 ,评价 
对 象 本 刁 存在 一 定 的 不 确定 性 ,各 种 复合 因子 相互 
影响 ,要 准确 客观 地 厘清 界限 是 比较 难 的 。 比 如 ， 
如 何 将 地 面 降水 量 中 人 工 增加 部 分 从 自然 降水 量 
中 识别 出 来 ,如 何 所 除 “ 假 效果 ”而 实现 精准 评价 以 
及 如 何 建 立 概念 模型 等 问题 *”。 区 域 历史 回归 统 
计 检 验 是 增 雨 效果 评价 的 基本 方法 ,但 往往 受 分 区 
分 型 方法 与 影响 区 、 对 比 区 设计 的 左右 ,很 显然 ,不 
同 对 比 区 将 会 生成 不 同 的 评价 效果 。 无 论 是 统计 
检验 ,物理 检验 和 数值 模拟 检验 ,在 对 比 区 选择 上 ， 


态 环境 可 持续 发 展 的 制约 因素 ,也 成 为 学 术 界 探究 
的 热门 话题 。 随 着 气候 变 暖 ,极端 天 气 影响 频繁 
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准确 性 是 很 难 把 握 的 。 要 科学 评价 增 雨 效果 ,精准 
分 区 尤为 重要 。 聚 类 检验 为 对 比 区 选取 提供 了 一 
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种 量化 方法 ,为 增 雨 效果 评价 的 定量 客观 提供 了 
科学 的 解读 依据 。 

多 年 来 气候 区 划 边 界 问题 因 受制 于 人 为 因素 
左右 ,其 地 理 分 类 方法 备 受 质疑 ,而 主 成 分 分 析 
(Principal component analysis , PCA ) 在 一 定 程度 上 规 
避 了 这 种 缺陷 ,但 在 荷载 值 上 依然 存在 临界 值 勘定 
问题 ,因而 在 一 定 意 义 上 尚 不 能 实现 科学 准确 的 数 
值 分 类 。 聚 类 分 析 又 称 群 分 析 ,是 分 类 学 中 依据 事 
件 的 某 些 相似 性 特征 进行 归 类 的 一 种 多 元 方法 ,在 
分 析 中 多 选用 距离 系数 、 相 似 系数 、 相 关系 数 、 相 关 
距离 系数 等 指标 以 实现 相似 性 度量 "”"。 尽 管 如 
此 ,很 多 学 者 依然 认为 该 方法 缺乏 统计 学 基础 。 乏 
枕 生 哺 提 出 的 聚 类 统计 检验 方法 (Cluster analysis 
with statistical test,CAST) 则 有 效 弥 补 了 分 区 分 类 过 
程 中 显著 性 检验 的 重大 缺陷 ,尤其 在 中 心 聚 类 评价 
上 具有 广泛 的 适用 价值 ”。 

根据 聚 类 统计 检验 分 析 的 基本 思路 ,本 文 对 祁 
连 山 及 周边 降水 空间 分 区 分 型 的 合理 性 进行 了 深 
入 探讨 ,人 研究 结果 将 有 利于 剖析 复杂 地 形 下 不 同 降 
水 区 域 的 空间 变化 特征 ,为 准确 把 握 祁连山 降水 变 
化 特征 与 分 布 规律 提供 科学 依据 。 


1 数据 与 方法 


1.1 数据 来 源 

本 文 数据 来 源 于 中 国 气象 局 国家 气象 信息 中 
心 提供 的 祁连山 以 及 周边 区 域 31 个 国家 基本 站 
1961 一 2020 年 5 一 9 月 降水 量 资料 ,研究 所 用 序列 经 
过 严格 质量 控制 。 如 遇 个 别 时 段 缺 测 , 则 按照 维尔 
TK (WILD ) 综 合 订正 法 取 邻 近 站 点 展开 序列 订正 ,在 
通过 适当 性 判 据 之 后 ,实现 该 站 点 缺失 值 的 插 补 55 。 
祁连山 及 周边 海拔 高 度 变 化 与 气象 站 点 分 布 见 网 
1。 祁 连 山 西高 东 低 ,最 大 海拔 高 度 差 2084 m, 为 了 
消除 地 理 因素 引起 的 平均 值 以 及 变 率 的 影响 ,在 分 
析 前 对 上 述 资料 进行 了 标准 化 处 理 ,形成 祁连山 降 
水 空间 场 数据 集 。 分 析 过 程 采 用 多 种 统计 方法 并 
借助 SPSS „Matlab 编程 以 及 ArcGIS 等 软件 实现 了 降 
水 特征 的 分 区 分 型 研究 。 
12 研究 方法 
1.2.1 主 成 分 分 析 ”特征 值 是 衡量 主 成 分 影响 度 的 
重要 指标 ,代表 了 引入 该 主 成 分 以 后 能 够 解释 源 数 
据 信息 量 的 多 少 ,一 般 以 特征 值 >1 作为 考量 主 成 分 
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图 1 祁连山 及 周边 海拔 高 度 变化 与 气象 站 点 分 布 
Fig. 1 Changes of altitude and spatial distribution of 
meteorological stations in Qilian Mountains and its 


surrounding areas 


重要 性 的 标准 。 方 差 贡 献 率 (或 累积 贡献 率 ) 用 于 
判断 主 成 分 的 份量 ,其 数值 越 大 ,表明 该 主 成 分 (或 
上 述 主 成 分 ) 包 含 的 源 数据 信息 量 就 越 大 。 和 荷载 值 
反映 了 主 成 分 与 各 个 变量 间 的 密切 程度 ,在 标准 化 
变量 下 ,任意 一 个 主 成 分 的 荷载 值 实质 上 就 是 主 成 
分 与 原 变量 的 相关 性 指标 “" ,也 就 是 说 ,在 一 定 意 
义 上 荷载 值 就 类 似 于 相关 系数 ,其 值 越 大 ,表明 该 
因子 对 当前 变量 的 影响 程度 越 大 。 

主 成 分 分 析 天 2 可 以 归纳 为 以 下 几 个 步 又: 
(1) 对 降水 量 数据 进行 标准 化 处 理 , 以 消除 地 域 因 
素 带 来 的 系统 性 误差 ;(2) 计算 标准 化 样本 的 相关 
系数 矩阵 ; (3) 按照 特征 值 和 方差 贡献 率 的 大 小 提 
取 因 子 数量 ,取得 阔 值 范围 内 的 主 分 量 ; (4) 对 各 个 
主 分 量 进行 分 析 ,考察 可 解释 性 ,如 有 必要 可 对 因 
子 进行 旋转 ,通过 荷载 值 以 寻求 更 好 的 解释 方式 。 
1.2.2 聚 类 检验 ” 聚 类 分 析 可 将 样本 划分 为 不 同类 
群 ,通过 评价 分 析 , 其 结果 具有 综合 性 .客观 性 和 科 
学 性 ,近年 来 在 各 个 领域 得 到 广泛 应 用 i。 气候 
分 型 区 划 理 论 研 究 已 经 证 明了 统计 聚 类 检验 与 主 
成 分 分 析 的 关联 性 ,并 认为 在 一 定 程 度 上 聚 类 检验 
是 主 成 分 应 用 于 分 区 分 型 的 理论 基础 ,并 通过 多 种 
方法 证 实 了 这 种 关联 的 有 效 性 ““”…。 聚 类 分 析 常 用 
TPP PS OR BE at OCH LB] AY SRK AB, PABA FT FE 
离 .曼哈顿 距离 . 欧 几 里 得 距离 SKIPS S. 
将 变量 间 这 种 距离 用 上 表示 ,Cy 与 相关 系数 具有 一 
致 性 ,同时 又 具备 正 态 分 布 性 ,通过 Xx 分布 函数 的 
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构造 ,可 以 实现 聚 类 过 程 的 显著 性 检验 。 
假设 某 气象 要 素 具有 时 间 序列 长 度 为 n. 测 站 2 结果 与 分 析 


数量 为 p 所 构成 的 变量 场 X, 已 经 证 明 , 测 站 闷 和 总之 
间 的 距离 可 以 用 相关 系数 (rj) 表 示 。 根 据 乏 枕 生 " 
理论 ,定义 自由 度 7-1 的 分 布 统计 量 : 


X=2n D0 cj 一 (1) 
i<j 

c= l-r; (2) 

w;=1/(1-r;) (3) 


6= Yoxe, / Qo; (4) 
i<j i<j 


式 中 : = fqq- 1) ,是 与 自由 度 相关 的 统计 量 ; 9 为 
同一 类 型 的 变量 个 数 , 对 于 气象 变量 场 X 来 说 , 9 就 
是 变量 场 中 同一 分 区 分 型 结果 中 的 测 站 数量 
(q<p) de i=l, 2, =, q- 1,j=1, 2, …, ;Gy 表示 距 
离 系 数 ;owi 表 示 距 离 权 重 ; 6 表示 同一 类 型 中 的 极 大 
似 然 估计 值 。 

根据 式 (1)~(4) 求 出 相应 站 点 妇 统 计量 ,在 给 
E fai BEa FW Lxi) A It HA, AE 
XY< 欠 -1), 则 表明 该 测 站 可 归属 于 同一 区 域 ， 
取 该 区 域 最 大 荷载 值 为 聚 类 中 心 , 所 有 符合 条 件 站 
点 即 为 同一 类 别 的 范围 。 
1.2.3 区 域 历 史 回归 检验 ”区 域 历史 回归 检验 方法 
是 基于 地 面 降水 量 的 效果 评价 方法 ,该 方法 立足 对 
比 区 .影响 区 地 面 降水 量 相关 分 析 ,采用 最 小 二 乘 
法 建立 一 元 线性 回归 方程 ,预测 出 影响 区 的 自然 降 
水 量 ,通过 自然 降水 量 估计 值 与 实测 降水 量 的 分 析 
比较 ,得 出 绝对 增 雨 量 和 相对 增 雨 率 。 该 方法 在 人 
工 增 雨 作业 效果 检验 技术 指南 中 已 有 介绍 , 详 见 文 
献 ,这 里 不 再 赣 述 。 


2.1 降水 量 分 布 基本 特征 

祁连山 及 周边 地 区 降水 量 分 布 与 祁连山 脉 走 
向 基本 一 致 。 以 5 一 9 月 为 例 , 东 段 降水 在 350 mm 
以 上 , 西 段 降水 只 有 100 mm 左右 ,表现 出 极 不 均匀 
性 。 降 水 量 年 度 分 布 不 均 , 以 1961 一 2020 年 为 例 ， 
全 域 多 年 平均 降水 量 343.7 mm, 降水 量 最 大 值 
441.8 mm(1967 年 ) ,降水 量 最 小 值 240.2 mm(1991 
年 ) ,降水 变 率 +30% ,高 值 区 多 在 祁连山 区 。 降 水 季 
节 性 特征 明显 ,夏季 集中 ,冬季 稀少 ,5 一 9 月 降水 量 
占 全 年 855% 以 上 。 降 水 地 域 分 布 极 不 均匀 ,以 2016 
年 为 例 ,张掖 降水 量 还 不 足 刚 察 的 1/45, 出 现 了 50a 
不 遇 的 气象 干旱 。 研 究 发 现 ,祁连山 及 周边 地 区 降 
水 量 与 地 理 位置 相 关联 ,降水 量 与 纬度 呈 显 著 负 相 
关 ,与 海拔 呈 显 著 正 相关 ,也 就 是 说 纬度 越 高 ,降水 
量 越 少 ,海拔 越 高 ,降水 量 越 多 ,这 种 规律 使 得 祁 连 
山 及 周边 地 区 降水 呈现 出 东 多 西 少 、 南 多 北 少 的 特 
征 , 最 大 降水 带 分 布 在 2500~3500 m 区 域 ,也 说 明 西 
南 暖 湿 气 流 与 夏季 风 对 该 区 域 降水 形成 有 着 极其 
重要 的 影响 。 
2.2 降水 量 的 聚 类 分 析 

建立 1961—2020 年 祁连山 及 周边 31 个 测 站 
5$ 一 9 月 地 面 降水 量 原始 矩阵 ,对 标准 化 后 的 相关 和 珑 
阵 做 主 成 分 分 析 ,然后 以 各 个 主 成 分 对 应 的 高 荷载 
值 作 为 聚 类 中 心 ,依据 式 (1)~(4) 给 出 一 定 信和 度 下 
的 分 区 检验 判 据 ,对 分 区 结果 进行 分 析 评 价 。 

首先 对 降水 量 相 关 和 矩阵 进行 主 成 分 分 析 ( 表 
1) ,依据 特征 值 以 及 方差 贡献 率 得 选 之 分 量 。 从 表 
1 可 以 看 出 ,旋转 前 后 7 个 主 分 量 其 累积 方差 贡献 率 
均 为 78.135% ,经 过 旋转 之 后 主 分 量 C:~C? 的 方差 贡 


表 1 祁连山 及 周边 降水 量 的 主 成 分 分 析 结 果 


Tab.1 Results of principal component analysis of precipitation in Qilian Mountains and its surrounding areas 


主 分 量 主 成 分 分 析 (PCA) 鲁 棒 主 成 分 分 析 (RPCA) 
特征 值 方差 贡献 率 /% 累积 贡献 率 /% 特征 值 方差 贡献 率 /% 累积 贡献 率 /% 
Ci 10.842 34.973 34.973 6.694 21.592 21.592 
C， 4.547 14.666 49.639 4.182 13.491 35.083 
G 2.942 9.490 59.129 3.824 12.334 47.417 
Cy 2.024 6.530 65.659 3.378 10.898 58.315 
Cs 1.512 4.877 70.536 2.257 7.282 65.597 
Co 1.341 4.326 74.862 1.959 6.318 71.915 
C 1.015 3.274 78.135 1.928 6.221 78.135 
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献 率 普遍 得 以 提高 ,说明 旋转 后 对 主 分 量 的 解释 更 
为 明确 ,这 7 个 主 分 量 可 以 很 好 的 表征 原始 场 大 部 
分 信息 和 主要 趋势 。 相 比 主 分量 C,, 其 余 主 分 量 的 
方差 贡献 率 都 比较 小 ,这 说 明 祁 连 山 及 周边 的 地 形 
影响 特征 很 明显 。 

以 主 成 分 分 析 得 到 的 高 荷载 值 作为 中 心 点 ,对 
答 载 值 > 0.6 各 个 测 站 进行 聚 类 评判 , 当 评 判 结果 通 
过 95% 显 著 性 水 平 检验 时 ,可 认为 分 类 结果 合理 ， 
评判 标准 与 结果 见 表 2。 通 过 上 聚 类 检验 ,人 研究 区 域 
构成 以 乐 部 、 海 受 、 野 牛 沟 武威、 高 台 、 临 泽 、 刚 察 
为 中 心 点 的 7 个 区 域 (Z1~Z7) ,一 个 区 域 对 应 一 种 类 
型 ,从 而 得 到 祁连山 及 周边 多 个 相对 独立 的 分 区 分 
型 区 域 。 
2.3 聚 类 结果 特征 分 析 

依据 聚 类 检验 分 析 结 果 ,对 各 个 降水 区 域 分 布 
特征 进行 研究 (图 2)。 结 果 显 示 : 第 一 区 域 Z1 位 于 
研究 区 域 东 南部 ,包括 乌 鞘 岭 OL ACE RAB OF 
安民 和 ,化 隆 、 循 化 8 个 站 点 , 占 研究 区 域 测 站 比例 
H 25.8% ,最 大 和 荷载 中 心 在 乐 都 (0.892 ) ,平均 荷载 
值 0.376 ,年 平均 降水 量 296.3 mm ,海拔 1800~3000 
m, 属于 人 研究 区 域 季风 气候 响应 区 ,是 祁连山 及 周边 
自然 生态 门户 以 及 农 牧 业 重点 发 展 地 区 。 第 二 区 
域 Z2 位 于 祁连山 中 南部 ,包括 门 源 , 海 受 EU .大 
通 .互助 .西宁 ` 涅 中 7 个 站 点 , 测 站 占 比 为 22.6% ， 
最 大 荷载 值 0.751 ,年 平均 降水 量 398.5 mm, 海拔 
2500~3000 m ,是 研究 区 域 草原 与 牧场 基地 ,也 是 祁 
连 山 区 重要 的 水 源 涵 养 地 。 第 三 区 域 Z3 分 布 在 沿 
祁连山 中 段 ,包括 山 丹 、 民 乐 . 肃 南 FE PAE 
祁 连 6 个 站 点 , 测 站 占 比 为 19.3% ,最 大 荷载 值 
0.786 ,年 平均 降水 量 285.0 mm ,海拔 1700~3400 m, 


是 祁连山 重要 生态 保护 与 植被 恢复 地 区 。 第 四 区 
域 74 分 布 在 祁连山 北 芳 东 端 即 石 羊 河流 域 , 测 站 占 
比 为 12.9% ,最 大 荷载 值 0.903 ,年 平均 降水 量 170.3 
mm ,海拔 1300~2000 m, 是 发 端 于 冷 龙 岭 北 侧 大 雪 
山 的 重要 绿洲 流域 ,属于 典型 的 灌溉 农业 。 第 五 区 
域 Z5 为 最 北端 琉 勒 河流 域 , 测 站 占 比 6.5%, 最 大 荷 
载 值 0.824, 年 平均 降水 量 仅 有 81.2 mm ,海拔 1400 
m, 是 源 于 祁连山 西 段 托 来 南山 与 跑 勒 南山 之 间 的 
绿洲 地 带 , 属 于 灌溉 农业 。 第 六 区 域 76 地 处 黑河 流 
域 ,发端 于 祁连山 北 项 中 段 冰川 地 带 , 测 站 占 比 为 
6.5% , 最 大 荷载 值 0.952 ,年 平均 降水 量 100.4 mm, 
海拔 1470 m, 属于 灌溉 农业 。 第 七 区 域 Z7 为 祁连山 
海拔 最 高 区 域 , 包 括 天 峻 、 刚 察 , 测 站 占 比 6.5% ,最 
大 和 荷载 值 0.814 ,年 平均 降水 量 341.3 mm ,海拔 3400 
m, 位 居 大 通 山 地 段 , 是 祁连山 多 条 河流 的 发 源 地 。 
综 上 可 以 看 出 ,分 区 分 型 结果 与 降水 量 地理 地 形 、 
海拔 高 度 显 著 相 关 ,各 个 分 区 相互 独立 ,差异 性 明 
显 , 同 一 区 域 特征 相近 ,具有 共性 ,从 而 构成 祁连山 
独特 的 降水 空间 布局 特征 。 

依据 聚 类 检验 方法 ,可 将 祁连山 及 周边 划分 为 
7 个 区 域 , 其 结论 与 研究 成 果 ”% 相 吻合 : 受 西风 带 、 
副热带 青藏 高 原 环流 的 共同 影响 ,特别 是 山地 地 
形 的 热力 和 动力 作用 ,祁连山 及 周边 容易 构成 气旋 
式 环流 ,副热带 高 压 明 显 增强 ,西南 暧 湿 气 流 强化 ， 
夏季 风 显 著 北 延 ,冷暖 空气 一 旦 相遇 , 极 易 在 祁 连 
山 一 带 形成 强 降 水 区 ,西北 气流 型 . 低 槽 东 移 型 AR 
高 西 低 型 均 是 构成 祁连山 强 降 水 的 有 利 天 和 气 形 
势 。 也 就 是 说 ,局 地 环流 .边界 层 结构 均 能 有 效 促 
使 祁连山 地 形 正 反馈 ,导致 降水 量 显著 增加 , 且 从 
东南 到 西北 逐渐 减少 ,从 而 构成 独特 的 降水 空间 变 


表 2 祁连山 及 周边 降水 量 聚 类 中 心 以 及 分 区 分 型 结果 


Tab. 2 Cluster center site and regional classification of precipitation in Qilian Mountains and its surrounding areas 


降水 区 域 中 心 站 最 大 荷载 值 站 点 数量 统计 量 临界 值 
Ži 乐 都 0.892 8 34.303 40.113 
72 海 受 0.751 7 7.607 31.410 
73 野牛 沟 0.786 6 2.054 23.685 
Z4 武威 0.903 4 0.501 11.070 
Z5 高 台 0.824 2 0.000 3.841 
76 临 泽 0.952 2 0.000 3.841 
Z7 刚 察 0.814 2 0.000 3.841 


注 :临界 值 是 指 在 给 定 信 度 0.05 下 具有 自由 度 w-1 随机 变量 x 的 检验 值 。 
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图 2 祁连山 及 周边 各 个 


区 域 荷载 值 的 空间 分 布 


Fig. 2 Spatial distributions of the loading values of each divided areas in Qilian Mountains and its surrounding areas 
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化 特征 。 
2.4 典型 流域 增 雨 效果 评价 
2.4.1 增 雨 效果 评价 党 景 ”统计 检验 的 评估 对 象 是 
地 面 降水 量 , 基 本 方法 通过 同一 时 期 对 比 区 自然 降 
水 量 与 影响 区 实测 降水 量 的 比较 ,得 出 定量 化 结 
果 ,这 个 结果 就 是 增 雨 效果 。 但 在 作业 影响 区 确定 
以 后 ,不 同 对 比 区 就 会 有 不 同 的 增 雨 效果 。 聚 类 检 
验 分 区 分 型 结果 ,可 以 作为 对 比 区 选择 的 基本 依 
据 , 换 一 句 话说 , 当 影响 区 和 对 比 区 位 于 同一 区 域 
时 ,统计 检验 才 有 了 可 信和 度 。 

选取 Z4, 即 石 羊 河流 域 作为 实证 区 域 。 该 流域 
人 口 密 集 、 水 资源 匮乏 、 生 态 环境 极其 脆弱 , 人工 增 
雨 已 成 为 缓解 旱情 与 生态 治理 的 重要 措施 。 现 已 
建成 人 工 影响 作业 点 71 个 ,其 中 火箭 增 雨 作业 点 45 
个 高 炮 增 雨 防 起 作业 点 14 个 固定 烟 炉 作业 点 8 
个 . 焰 弹 作业 点 4 个 ,覆盖 面积 6000 km 以上。 该 流 
域 人 工 影响 作业 始 于 1992 年 ,截至 目前 ,累积 作业 
6000 Arve ,发 射 高 炮 25000 余 发 .火箭 16000 余 枚 、 
焰 弹 20000 余 枚 ,具备 对 比 区 与 影响 区 遂 选 条 件 。 

首先 ,Z4 所 在 区 域 500~700 hPa 常年 盛行 西北 
风向 ,只 要 对 比 区 位 于 影响 区 上 风 方 或 者 侧 风 方 ， 
就 可 以 使 对 比 区 不 受 催化 作业 影响 而 保持 相对 的 
独立 性 ;其 次 ,Z4 所 在 区 域 地 处 大 陆 腹地 ,气候 干 
FA 降水 少 、 燕 发 量 大 ,地 形 地 貌 .气候 背景 、 植 被 特 
征 、 区 域 面积 .天 气 系统 高 度 一 尾 , 地 面 降水 量 显著 
相关 。 石 羊 河流 域 地 理 地 形 以 及 增 雨 作业 点 布局 
见 图 3。 
2.4.2 增 征 效果 评价 和 检验 ” 石 羊 河流 域 地 面 作业 
集中 在 5 一 9 月 ,选取 民 勤 为 对 比 区 , 凉 州 为 影响 区 ， 
即 对 比 区 位 于 影响 区 侧 风 方 。 将 1961—1991 年 作 
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图 3 石 羊 河流 域 地 理 分 布 以 及 对 比 区 与 影响 区 布局 
Fig. 3 Geographical distribution of Shiyang River Basin and 


the layout of contrast area and affected area 


为 历史 降水 量 序列 时 段 ,1992 一 2020 年 作为 作业 影 
响 序列 时 段 ,2 个 时 段 样本 量 基本 相当 。 在 应 用 区 
域 历史 回归 统计 检验 方法 之 前 ,需要 确定 统计 变量 
是 否 满足 正 态 分 布 ,为 了 使 回归 模型 具有 可 比 性 ， 
要 求 对 比 区 和 影响 区 在 2 个 时 段 同 时 具备 正 态 性 ， 
即 通 过 Kolmogorov-Smimov (KS) 检 验 1。 检 验 与 分 
析 结 果 见 表 3。 

以 对 比 区 降水 量 为 自 变 量 ,影响 区 降水 量 为 因 
变量 ,运用 地 面 降水 量 区 域 历 史 回 归 统 计 检验 方法 
对 该 时 段 降水 量 进行 分 析 。 对 应 时 段 回 归 方 程 均 
通过 95% 的 显著 性 水 平 检验 ,表明 拟 合 方程 可 以 代 
表 该 时 段 降 水 量 的 基本 特征 ( 表 3)。 将 对 比 区 降水 
量 代 入 回归 方程 ,得 到 影响 区 降水 量 的 预 估 值 ,在 


表 3 石 羊 河流 域 不 同时 段 地 面 作 业 效 果 检 验 的 特征 值 


Tab. 3 Eigenvalues of ground operation effect test at different periods in Shiyang River Basin 


Z(Sig.) 


月 份 相关 系数 (7) 回归 方程 7 ”绝对 增 雨量 /mm ”相对 增 雨 率 /% 
对 比 区 影响 区 
5 1.015(0.254 0.935(0.346) 0.833" y=5.929+0.950x 8.110 -0.19 -1.04 
6 1.242 (0.092 1.257(0.085) 0.604 y=10.411+0.972x 4.084 -1.56 -6.26 
7 0.864(0.444 0.905(0.386 ) 0.490" y=18.367+0.385x 3.029 6.30 21.86 
8 0.948 (0.330 1.051(0.220) 0.680" y=12.292+0.832x 4.998 5.44 16.11 
9 1.261 (0.083 0.877(0.426 ) 0.677 y= 9.889+0.849x 4.948 1.26 4.40 
5 一 9 0.502(0.963 0.564(0.908 ) 0.665" y=54.432+0.800x 4.798 8.91 6.51 


LE : Kolmogorov-Smimov(KS ) 非 参数 检验 中 Z 为 检验 统计 量 , 当 双 尾 渐进 概率 Sig.>0.05 时 ,说 明 样本 具有 正 态 性 ;* 表 示 相 关系 数 通过 95% 显 
著 性 水 平 检验 ;7 为 线性 回归 方程 统计 检验 值 ™, 当 70.367 时 ,认为 回归 关系 通过 95% 显 著 性 水 平 检验 ;绝对 增 雨量 = 影响 区 降水 量 一 对比 区 
降水 量 ; 相 对 增 雨 率 =( 影 响 区 降水 量 /对 比 区 降水 量 - 1)x100%。 
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此 基础 上 对 不 同时 段 作 业 效果 展开 分 析 。 

1992—2020 年 5—9 月 绝对 增 雨 量 为 8.91 mm, 
相对 增 雨 率 为 6.51% ,说 明 5 一 9 月 作业 效果 明显 ， 
具有 很 好 的 增 雨 潜力 。 但 在 不 同月 份 作业 效果 却 
表现 出 明显 的 差异 性 ( 表 3)。5 一 6 月 为 负 效 果 ,7 一 
9 月 为 正 效果 ,其 中 7 月 最 高 (6.30 mm, 21.86% ) ,8 
月 次 之 (5.44 mm, 16.11%) ,这 一 结果 与 石 羊 河流 域 
气候 背景 以 及 天 气 系 统 演变 规律 密切 关联 。5 一 6 
FATE AB AR AR SW), ,是 典型 的 “ 卡 脖子 ”时段 , 从 多 年 
降水 量 分 布 上 看 ,该 时 段 降水 量 不 足 全 年 的 20% ， 
水 汽 条 件 差 . 增 雨 效率 不 高 .作业 效果 弱 ,实际 业务 
中 应 该 慎重 选择 作业 时 机 , 尽 可 能 降低 负 效 果 的 出 
现 几 率 。 而 7 一 9 月 正 是 高 原 降水 云 系 活跃 阶段 , 尤 
其 是 7 一 8 月 , 阵 性 降水 多 .水 汽 条 件 好 、 增 十 效率 
高 ,作业 效果 明显 ,可 以 通过 增加 作业 点 次 以 及 作 
业 量 ,显著 增加 降水 量 , 有 效 提升 作业 效果 。 


3 讨论 


作业 效果 评价 一 直 是 制约 人 工 影响 天 气 发 展 
的 技术 性 难题 。 祁 连 山 地 形 复杂 、 范 围 广大 ,我 们 
前 期 已 经 对 降水 分 布 的 区 域 性 特征 进行 了 一 些 研 
究 , 在 利用 统计 检验 评价 增 雨 效果 时 ,发 现在 对 比 
区 的 选择 上 存在 诸多 问题 ,导致 检验 效果 出 现 很 多 
不 确定 性 结果 。 众 所 周知 ,对 比 区 与 影响 区 降水 量 
比 对 是 统计 检验 的 基本 思路 ,大 量 事实 证 明 ,在 作 
业 影 响 区 确定 以 后 ,不 同 对 比 区 就 会 有 不 同 的 增 雨 
效果 。 遗 憾 的 是 ,一 直 以 来 ,关于 对 比 区 选择 并 未 
形成 客观 准确 的 方法 ,也 没有 可 以 依靠 的 理论 依 
据 ,地 面 作 业 效 果 评价 与 人 工 增 雨 服务 陷 人 困境 。 

聚 类 检验 分 区 分 型 结果 ,可 以 作为 对 比 区 选择 
的 基本 依据 。 只 有 当 影 响 区 和 对 比 区 位 于 同一 区 
域 时 ,统计 检验 才 具 有 可 靠 的 前 提 条 件 。 本 文采 用 
主 成 分 分 析 方法 对 研究 区 域 降水 场 进 行 分 区 分 型 ， 
再 通过 聚 类 检验 进行 结果 验证 ,很 大 程度 上 保证 了 
分 区 的 准确 性 。 将 分 区 分 型 结果 应 用 于 对 比 区 与 
影响 区 的 选择 ,解决 了 在 增 雨 效果 评价 上 区 域 选 择 
的 盲目 性 问题 ,这 也 是 制约 地 面 增 十 效果 评价 的 关 
键 环 节 。 这 种 方法 通常 用 于 气候 分 型 区 划 ,但 很 少 
用 于 区 域 降水 研究 ,在 人 工 增 雨 效果 检验 对 比 区 选 
择 上 则 更 为 少见 。 祁 连 山 以 及 周边 区 域 地 形 复 杂 ， 
降水 分 布 极为 不 均 ,建立 科学 准确 的 分 区 方法 ,对 


把 握 该 区 域 降水 特征 具有 极为 重要 的 意义 。 


4 结论 

基于 主 成 分 分 析 与 聚 类 分 析 的 一 致 性 ,本 文 对 
分 区 分 型 方法 进行 了 深入 研究 ,并 验证 了 这 种 结 
的 合理 性 。 通 过 分 析 ,将 祁连山 及 周边 降水 分 布 划 
分 为 7 个 区 域 , 在 对 上 述 区 域 逐 一 评价 的 基础 上 , 选 
择 石 羊 河流 域 作 为 实例 。 设 计 对 比 区 .影响 区 , 采 
用 区 域 历 史 回 归 统 计 方 法 对 不 同时 段 降水 量 展开 
增 雨 效果 评价 。 结 果 显 示 :1992 一 2020 年 5 一 9 月 绝 
对 增 雨 量 8.91 mm, HAXHIN YK 6.51% ,其 中 7 月 最 
高 (6.30 mm, 21.86%) ,8 月 次 之 (5.44 mm, 16.11%), 
说 明 人 工 增 雨 效果 在 不 同时 段 具有 差异 性 ,适时 作 
业 、 科 学 作业 ,才能 提高 综合 效果 ,最 大 程度 提升 抗 
FEK .生态 治理 与 修复 的 能 力 和 水 平 。 
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Cluster analysis with statistical test of precipitation distribution in Qilian 
Mountains and its surrounding area and evaluation of artificial precipitation 


enhancement in typical watershed 


GUO Xiaoqin'”, LI Guangming’, SUN Zhanfeng’, WANG Xingtao* 
(1. Key Laboratory for Cloud Physical Environment of China Meteorological Administration, Beijing 100081, China; 
2. Wuwei Meteorological Bureau, Wuwei 733000, Gansu, China) 


Abstract: The topography of Qilian Mountains and its surrounding areas in northwest China is complex, and the 
precipitation it receives is extremely irregular in time and space. The ability to characterize the distribution of 
rainfall and snowfall has great practical significance for understanding changes in precipitation patterns. Using 
data from 31 base stations covering the months of May-September, 1961—2020, precipitation was analyzed by 
principal component analysis (PCA). The significance of the results was then tested by cluster analysis with 
statistical test (CAST). Finally, a geographic zonation was applied to evaluate the effects of artificial precipitation 
enhancement. The results show that: (1) Seven geographic areas (Z1- Z7) can satisfactorily represent the 
characteristics of precipitation distribution in the study area (the contribution rate of cumulative variance exceeds 
78%). (2) The zones take Ledu, Haiyan, Yeniugou, Wuwei, Gaotai, Linze, and Gangcha stations as the central 
points respectively, and their zonal characteristics significantly reflect precipitation, terrain, and altitude. (3) 
Effects of artificial precipitation enhancement were evaluated by statistical regression of regional historical data. 
For the months of May-September during 1992—2020, the benefit was most pronounced in the Shiyang River 
Basin, where the absolute enhancement was 8.91 mm and the relative precipitation enhancement rate was 6.51% 
with the highest in July (6.30 mm, 21.86% ), followed by August (5.44 mm, 16.11% ). Assessing artificial 
enhancement of precipitation based on ground data relies heavily on the selection of a suitable comparison area, 
so results can change significantly with different comparison choices. The method of cluster analysis with 
statistical testing is not only helpful in zoning and classification of precipitation in complex terrain, but it also 
provides an objective basis for the selection of comparison areas. 

Key words: cluster analysis with statistical test (CAST); principal component analysis (PCA); division and clas- 
sification; evaluation of precipitation enhancement effects 


